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ABSTRACT 


The thesis covers the mechanical, thermodynamics, and aerodynamics 
of designing a free vortex axial flow turbine with the aid of a digital 
computer, It uses the Improved Ainley-Mathieson Performance Estimation 
Method to determine the aerodynamic pressure losses so off-design and 
design point catculations can be conducted to produce performance curves. 
It also gives physical and thermodynamic properties of the turbine. 
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Introduction 1 


fhe purpose of this thesis is to develop a computer program 
which can be used in the Preliminary Design of a Free Vortex Axial 
Flow Turbine. fhe design plan consists of three main parts. The 
first will be cycle calculations to determine stage characteris- 
tics at each station such as angles, density, temperature, flow areas, 
etc. The properties from this phase will ve used in parts 2 and 3. 
The second part is to determine mechanical characteristics Such as 
size and weights of turbine parts and stress and vibrational char- 
acteristics. The final phase is to determine turbine characteris- 
tics at design and off design parameters such as speed, pressure ratio, 
etc. This last phase is optional since compressible flow tables and 
gas properties are required and many turbine designers may be more in- 
terested in the size and weights than off desien characteristics in 
the preliminary design phase. These three parts are covered in cha- 
pters one througn three, chapter four will cover the operation of the 
computer program developed in this text. Appendix one is a list of 
symbols used in the text, appendix two is a user's guide to the com- 
puter program and has an example turbine design which is compared 
with a design proposed by reference 20. Appendix three is the com- 
puter program. 

The computer program is set up So it has many options of in- 


put variaoles so that as the designer becomes more knowledgeable 
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in his design study he can specify more characteristics and have 
more control over the design output. There are checks in the program 
to limit poor design properties such as excessive swirl in Output, 
high tip to hub ratios, negative reactions, etc. The computer pro- 
2ram can be used to design a constant hub, mean, or tip turbine. The 
design does not consider the effects of cooling in turbines, Since 
time was not available to incorporate this into the computer program. 
This limitation mostly affects aircraft gas turbine and aircraft de- 
rivative marine gas turbines since these engines operäte at hish tur- 
inlet temperatures and cooling is required to allow for a high mean 
time between failures. 

There has been no effort made to prove the relations used in the 
design calculations. These can be found in basic physics, thermody- 
namics, and fluid mechanics books, 

The notation as to signs of angles in the velocity diagrams will 
be the same as used in chapter 4 of Gas Turbine Handbook which is re- 
ference (3). An illustration is enclosed in chapter | of the sign 
convention and notation used in this text. A deviation from this no- 


tation will be used in section (3.2). 








Chapter | 9 


Cycle Calculations 


I The basic design of an axial flow turbine starts with the 
cycle calculations to determine the thermodynamic properties thro- 
ughout the turbine. The basis of these calculations are from known 
inlet conditions and assumed design parameters. The computer pro- 
gram is set up so that a minimum number of these parameters must be 
specified and leaves the option to the designer to more fully spec- 
ify the design parameters, if he so desires. This chapter intends 

to show how these basic calculations are carried out and the methods 
employed to determine the physical properties such as flow area from 
these design parameters. Figures (1-1), (1-2), and (1-3) show the 
sign convention and notation that will be used in this chapter. The 
following are some of the design parameters nceded for the cycle cal- 
culations; inlet temperature (^R) Toi, turbine inlet pressure (psi) 
Poi, omega (rpm)/2 , specific heat (btu/lbm °R) Cp, mass flow 
(1bm./sec)m, ratio of specific heat SX. tip radius (ft)rr,tip to 
hub ratio (rt/rh ), desired efficiency/l, loading coefficient P, etC. 
Appendix two (2) will specify which parameters must be known and 
which ones are optional. 

122 Determination of overall turbine parameters fron inlet 
conditions And assumed design parameters. 


(a) (152590) 
Power (hp) 550 
Ze - Power (hp) 550 


mJ ما‎ (BE 


Toe = Toi . 2۶ Cr) 
Cp 
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Fig. 1-3 T-S diagram for a Reaction stave 


Rauyeroplic Efficiency 11 


Np 8 In(Toe/Toi) | (ООО т) 
In | (Toe/Toi-l) + 1 
Ü 

Poe = Poi (Toe/Toi)’PpTy=1) (1.2-4) 


(b) If the pressure ratio were specified insteao of ootput power 


the following calculations will be performed. 


Polytropic Efficiency X.1 

Cf Ln(1 -JE(1-(Poe/Poi 7)) — 

Ink Poe/ Poi aS 

Y 
Toe = Toi (Poe/Poi)"fp(X- 1) (1.2-6) 
/Aho = Cp(Toi - Toe) EE 
Power out = Anodes __ (1.2-3) 

550 
(c) 

رو (rad/sec)‏ جا سوہ 
n D CAM (1.2-10)‏ 
)1.2-11( 2 طض TE (CT‏ 


Un S r aN ft/sec) (1.212) 





Euler's Equation 7۲ 


| UhVx (tai? - tan) 
20 J 
coJ 人 ho 
Y E eX stage (1.2-14) 
2 
Uh 


The value of Wis taken at the hub or the last stage because this 
is the most critical in the design of a turbine in reference to load- 
ing. Normally the value of P is from (1.5 - 2.9). The lower value is 
for lightly loaded stages while the higher is for highly loaded stages. 
To have no exit swirl Y must be less than or equal to 2,0. 


Ye2.0 ==> (no exit swirl) 


for (X= Ú 5 
h R= Pr diss 955 


where 
ba > ha ratio of change in static enthalpy 
(R= across the rotor to change in (1.2-15a) 
hoz- hos 


stagnation enthalpy across the stage. 
With the value of J speci fied, the number of stages required to 


develop the power out desired can be determined. 


2 
Aho, taze 7 pus O 


goJ 


Determining number of stages 
nor 


NOs taqe 


N = integer 


If N were specified instead of Y, the change in enthalpy pet 


stage would be determined. 
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13 
This program uses a constant change in enthalpy for each stage. 


1f 70 and X h = 0 the radius of the hub can be determined from 


equation (1.2-17). 


(1.2-17) 





Another parameter which will be used is D = Vx/Uy (flow coeffic- 
ient). Where Vx is the axial velocity, which is constant in a free 
vortex turbine. This parameter ranges from 0.4 - 1.3, the higher 
values are used for aircraft turbines which use high axial velocity 
for propulsion. [In marine and industrial turbines a low value is 
desired because the energy of the exhaust gases are not recoverable 
in this form, unless it is used to drive another turbine. 

Another parameter that may be specified is the critical Mach 


number at the exit, 


VX 


2 2 
En m 
where R is the gas constant 
R = Cp(X-1/ HJ روحم‎ 


The Vx is computed exactly 


Vx = Mo. Keraosonme | (1.2-13a) 
Crit 
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If the exit Mach number is specified the turbine exit area is 


varied until the required.area is obtained to give the specified ax- 


ial velocity Vx. If this number is not specified but the exit area 


is specified the value of Vx is varied until the required area is 


achieved. 


From the continuity equation for gases 


Vx = m (1.2-20) 
— 
2/7 POT 
r 
where h 
£o "ET шы a) 
RT(r) 
2 
T(r) = To - m (1.2-22) 
2 CpgoJ 
Y 
P(r) = Po/(To/T(r)) 8-! (1,2-23) 
V(r) (1225) 
r) = -一 一 一 一 2g 
cos (Y (r)) 


where r varies 


equal to zero 


hu be = 


between the hub and tip. If the exit angle is not 


tan Ha -(Rhe Дуза] (1.2-25) 


and for free vortex 


rVg = rVx(tan(Y ) = constant (1.2-26) 


then 


Kr.) 


- a O (12527) 
Ee 


where re is taken at hub mean, and tip. This calculation is canriad 


out in function progran ANGLE. Thus the exit angle has been deter- 


mined at the hub, mean, and tip radii. 
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The pressure, density, and temperature calculations are per- 
formed in subroutine TPL. The integral calculation is performed in 
subroutine DENR which uses Simpson's Rule to compute the axial vel- 
ocity Vx. with the absolute angles(X(re) determined the relative. 
angles (rg) can then be determined by equation (1.2-28). This is 


carried out in function program ANGLES. 


Sf U 
Bro) = tan entsat - =. (1.2-28) 
VX 
From equation (1.2-15a) 
N لد للك‎ _ \хсаафы)_ (1.2-29) 
一 -二 一 - 一 一 
20020 Wr, 


which is computed in function prozran REACTI. With the exit cond- 
itions specified, the properties at station | and 2 can then be 
determined for each stage. 

The following properties are now known at the exit of the turbine; 
Vx, Toe, Poe, fra), PG), TG, (Ge, OE апа с. ‚ where 


r 


= r r ER 
e hub? “mean? tip 


13 Determination of properties at station 1 for stage n where 


۲92 1 « « < N. 


103 + Ahostase _ CIO) 


To} = 
Cp 
Bot oa € Po34/(Toagn/ 0ج‎ 29080 (1.3-2) 


۶۰٢٠٠1 lst stas Xin = 0.0 where n = 1 
TE hot ON, air. = Cadre! where re is the radius vhich remains 


constant throush out the turbine, ie(constant hub, tip, or mean). 
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16 
tan asao - ке (1.3-3) 


ED 





۷× 

2 
Tınlı c) = Tojn - Vif rc) ano) 

220JCp 

TE VX 
一 (15555) 
_ cos (jnre) 
e 
Pin(r.) = Posn/(Toyn/Tın) El (13-6) 
"ту CREARE c 0 
لك‎ RT;n(rc) 


which is identical to equations in section for exit of the turbine. 


A trial blade height (h) at station In will be established. 


o. (1.3-8)‏ ٭ مہہ 


2 xr e (Cre) 
Let delta һ = h/10.- Now depending on whether the turbine is constant 


hub, mean, or tip, radius are established at the other radii. 


Constant hub; 


Кп паб)” Te (1.3-9) 
r ¡¡ (mean) کے‎ h/2 (1.3-10) 
rin(tip)=- rc + h GERS 


Similarly for constant mean or constant tip. 

Of can then be determined at all radii for station l and temo- 
erature, pressure, and density can be calculated at all radii in the 
same manner as was carried out at the exit of.the turbine. Then 


Vx prime can be calculated in the same manner as Vx was calculated. 
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Vx prime = CIS 
Dm 
When Vx prime = Vx the proper radii has been calculated at station 

L. The method employed continuously decreases delta h until conver- 
zence is accomplished. This is the same method which was used at the 
exit of the turbine and will be used in station 2 calculations. The 
properties ot the gas is now determined as per exit of the turbine for 
station l. If tnis isn't the first stace the properties at station 1 
are the same as station 3 of the preceedinz staze, 
1.4 Determination of properties at station 2 for stage n where n = 


1 ° ۰ °” o N 


Toan = Tola (1.4-1) 

onlo) = tan [goJAhos د‎ tan(O@n(r.)) (1.4-2) 
۷× WF « 

EE vx/cos(Kan(r.)) ANS 
2 

Tan = Toan - V(rc)/220JCp (1.5-4) 

Pan sei P4n( T2n/ 4n) "1p CY-1) (Lab) 


which is an approximation from figurc (1-3), where the actual pressure 
is determined from pressure loss data as will be covered is Chapter 3. 


Pom = a! (1.4-0) 


fe = E: PN (1.4-7) 
RT2n(r) 
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These calculations are accomplished in subroutine TPD2 to 18 


determine temperature, pressure,and density at station 2. The calcu- 
lations are carried out in the same manner as station 1 calculations , 
to determine the gas properties and radii st station 2. The calculations 
are then carried out at station l of the next stage(ie; n = n-l until 
properties at all stations and positions have been computed for all 
stazes.) 

8 Having solved for gas properties at all stations and positions 
and having flow areas and gas angles the relative stagnation temperature, 


and the relative Mach numbers for stators and rocors can be calculated, 


we = 
Tree I, + (1.5-1) 
2g0CpJ 
L 
2 
Nnozzbe ‚per 8 (Toin/Tın) - DA (1.5-2) 


5 
anoon/ Ten) - 1) ge (1.5-3) 
^ 


D 


Mrotor = (Tozn,.1/T2n) - Y = (1.5-4) 
inlet X- 
D 
2 2 
ےہ7۸‎ Kran Ten) - 1 (195-5) 


These calculations are carried out at each stage. The followin? 
parameters have been determined at all stations of the turbine; temn- 
erature, pressure, density, flow area, absolute and relative anales, 
Reaction, relative Mach numbers, and the axial velocity Vx. These 
properties will be used to determine the mechanical aspects of the 


Turbine design and be used in performance estimation. 
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Chapter 2 19 
Mechanical Design of an Axial Flow Turbine 
2-1 In this chapter the mechanics of the program and assumptions 
made to design nozzle and rotor blades, and the turbine disc will be 
discussed. 
202 Having the gas properties and flow angles throughout the tur- 
vine the blade angles can be calculated for the turbine. The criter- 
ia used in establishing the blade angles is the loss data provided by 
figure (2-2) and (2-3) based upon Ainley Mathieson method of perform- 
ance estimation. As can be seen for reaction blades in figure (2-2) 
the incidence angles can vary from -15 to + 15 without much of an in- 
crease in profile loss Yp. In actual construction this large varia- 
tion can allow for a decrease in twist in the blade from hub to tip 
normally required in free vortex blades. For the computer progran, 
selection of blade angle, at the inlet of the blade is for zero incid- 
ence at design point. 

Figure (2-3) shows how blade exit angle varies with gas angle 
there is also a correction for Mach number. For Mach numbers less 
than 0.5 the values given by figure (2-3) will be used, for values 
greater than Mach 1 the blade angles will be equal to gas angles 
while for values between Mach 0.5-1.0 a linear variation will be as- 
sumed as is shown in figure (2-4). The calculation of these blade 
angles are carried out in subroutine Blade, which also will be used 
to determine gas angles when blade angles are known for off-design 


calculations in chapter 3. There is an additional correction for blades 
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Figure 2-1. Conventional Blade 
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Figure 2-2. Variation in profile 
loss with incidence angle( From 
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Figure 2-4. Angle variation 
with mach number( From D.G.Ainlev 
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21 
with a curved back trailing edge in reference (2) but this will not 


be applied in the computer program. 

2,3 The computer program can now determine the ideal pitch to chord 
ratio P/C, blade chord (C), aspect ratio H/C, bending and centrifugal 
stresses, number of blades, lst bending frequency, and weight of blades, 
As will be shown in this section all these parameters are interrelated 
and they all help determine each other. 


2.3a Centrifugal Stress Calculations 





Figure 2-5. Rotating Turbine blade 


dE e dmr u) = dmr (2 1h0,/60)* (2:321) 
AGE) = Ab f(r) (2.3-2) 
dm = Ab f(r)dr -É- (2.3-3) 


dF = (2m 60). A(r)rdr 
¡30 
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profile area ratio (Fron Emmert, Current Desian 
Bractices.for Gas Tucbine Elewents.(9)) 


2 
ES THOME дд... (2.3-5) 
$ 1000 


where; Ec = average tensile stress(psi) 
Tf = taper factor(from fizure (2-6)) 
5 specific mass of the blade 
A = annular flow area of blade ring (sq. in.) 
Ad. = shaft speed (rpm) 
For the computer program Tf will be based on linear variation with 
tip to hub area ratio (arcara) as an input, normally this will be 


between J.25 - 0.3233 
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[f a linear taper is assumed equation (2,3-1) can he integrated 


directly assuning; 


Alr) = Ab(1 - Xr) (2, узб) 
С1р 
دہ‎ fa =X r)rdr (2.3-7) 
Chub 


2172 Y X - rfub ) X,z3 3 
(ше go وم‎ ee ae = гер? 


2 
E ہ2114‎ (UE: Be бш) (X 2 3 
= > Ab -@ Zl: 2 +3 hub” Bio 


/ (2.3-8) 
2۵00 can be evaluated 
Atip = DN area ratio (Ab) 
X = l - area ratio (2.3-9) 
"tip 
229b Determination of ideal pitch to chord ratio. This value for 


P/C ratio is taken from figure (2-7) which is a correlation of Ainley's 
optimum spacing data. For stators the inlet gas angle to enter figure 
(2-7) is - OC, and the relative efflux angle is@, For rotors f is 


the inlet gas angle while -各 ;is the relative efflux angle. This cal- 
culation is done in subroutine STOCRA and the input anzles are taken 
at the mean radius. 

2. 3C Determination of non-dimensional section modulus(sM) and non- 
dimensional area at the base of the blade (Area/C“). An input to eet 
these values is the turning angle at the base of the blade. ¡e(hub 
for rotor and tip for stator). Figure (2-8) taken from reference(7) 


illustrates how the non dimensional base area is obtained. The non- 


dimensional section modulus is obtained from figure. (2-9) which is 


= 


50 2 un 


341332 
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Figure 2-7. Optimum spacing of turbine 
blades. (From a correlation of Ainley's 


loss data (2)) 


an unpublished relationship provided by Ainley for 
4 modification was made to this relationship based 
performed by Professor A.D. Carmichael to make the 
larger to conform to modern turbine blades. The SM 
2.35 for use in the computer program, 


2.3d 


Determining bending stresses on turbine blades. 


use in reference (6). 
upon calculations 
section modulus 


was multiplied by 


As can be seen 


in figure (2-10) and equation (2.3-10) a detailed calculation would be 


needed to determine the gas bending stresses on a turbine blade but 


because angle g is small and Mw is by far the greater bending moment an 


approximation proposed in reference (6) is used in 


design. 
D Vx (tanoQ- tamX3 Jh 
(@ gb) ax ای تام‎ 
2 ۵ ] 2 


(rotor) 


the preliminary 


(2.3-11) 


Where the angles are taken at mean radius and the following are the 


definition of the parameters; 
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Figure 2-3, Approximate Areas of Turbine 3lades vs 
Turning angle. (From paper by R.E.Dundas in Gas 
Turbine Handbook(7)) 
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Figure 2-9. Section Modulus as a function of blade 
turning angle. (From unpublished paper by Ainley 
for use in reference (6)) 
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X = S. (Ma cos - Mw sin ) 
Ogb- i. $ 7 
Y E 
+L (W cos + Ma sing ) 
[xx $ 7 
(2.3-10) 
Y 
Axial 
Š Direction 
Figure 2-10. 3ending moments acting 
on a turbine blade. 

Z = 3M C Section Modulus 
SM = non dimensional section modulus 
a æ= number of blades 

en where rg, = mean radius‏ يات 

Pitch 

e m Vx (tan Ko- tan oq) h 

o = - < 

ner 2 30 2 My SM c? sa 

h Vx (tan - tano) һ Р 
^" 0۲و‎ 5 C2 C 


At this point the blade chord has not been determined so the bending 
stress cannot be calculated, but equation (2.3-1la) will be coupled 
with a natural frequency equation and a cyclic loading equation to 


calculate the blade chord. 
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Figure 2-11 Effect of taper a 


frequency (Fror Emmert, Current Design Practices 
for Gas Turbine Power Elements (9)) 
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Figure 2-12. Variation of rotation coefficient with 
blade height ratio.(From Emmert, Current Desien 
Practices for Gas Turbine Power Elements.(9)) 
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250 Determination of natural bendiny frequency of blades from 23 
reference (9). The following expression is used in the computer 


program. 





k ME C (2.3-12a)‏ رت را رح 
a‏ 


f, is standing frequency (cps) 
Cf frequency correction factor taken from figure (2-11) 
k minimum radius of gyration of base profile(in.) 
h blade height (in) 
تد‎ ۱ 6 
E modulus of elasticity (psi) (29 - 30) x 10 psi 


ek specific weight of blade material (1b/in?) 


eo E E (2.3-12b) 
E S AS 
= 11.0 EE | tE assuming I «Zt 
3 Ag 


A cross sectional area of hase of blade (in?) 


t maximum thickness of blade (in) 


= (t/C) C 
“ع رك یں وو‎ 
t. «11.0 ££ LSM CXt/OE C 
7 h2 | (a/c?) 
= 11.0 Cf 20 0 5 chord (2.3-12c) 
h (А/С?) 


The running frequency f given by equation (2.3-13) is a function 


of the natural frequency and the spced of rotation of the blade. The 
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factor 3 is given by figure (2-12), A plot of this frequency as a 


function of operating speed is illustrated in figure (2-13) which 

is called a Cambell Diagram. Normally for good design practice the 
running frequency should avoid the lat four haramnics of the operating 
speed in its operating range. Far the design of the blade chord iin 
the computer program the ¡effect of the: rotatimnal speed as indicated 
in equation (2.3-13) will be neglected but á more stringent critenia 
using the harmonics of. the operating speed to estimate the ivibration- 


al stress in the bLade will be used. 


E + 717 (213-13)‏ سو 
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Figure 2-13. Cambell diagram, plot of frequency 
of vibrations as a function of operating speed. 
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CE Determination of: turbine’ blade óharacteristics básed on : 


cyclic loading requirements by use of a Goodman Diagram of steady 

and vibrational stresses. In reference (7) a use of the Goodman 
diagram is used to determine the chord of a compressor blade. This 
can also be applied to turbine blades. A Goodman diagram is a plot 
of vibrational stresses as a function of steady stresses. The vibra- 


tional stress is taken at 107 cycles at the operating temperature, 
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Figure 2-14, Goodman Diagram of Vibrational 
stresses as a function of steady stresses. 
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This plot can be approximated as a linear function. 
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Figure 2-15. Linearlized Goodman Diagram 


In reference (7) the vibrational stress is proposed as a function 


of the bending stress as shown in equation (2.3-14) 
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Figure 2-16. Amplification Factor. (From 
Paper of Trumpler & Owens (19)) 


where Amplification factor is shown in figure (2-16) as a function 


of the harmonic of operating speed. The 1.3 is a stress concentration 


factor specifiod in reference (7). From personal conversation with 
Professor A, Douglas Carmichael at M.I.T. the amplification factor 


can be approximated by 41/n where n is the harmonic of the operating 


speed. 
ff Y 60 (2.3-15) 
ee vega) 
From figure (2-15) 
Qvib = Cl - a steady (2.316) 


where ei steady = d centri fusal + CÀ sas bending 


and Cl is.the upper limit for vibrational stresses and C2 is the 
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upper limit for steady stresses, 
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ca ab 1.39760, XC (ger) 
í C2 C2 
Cl 1.3 (42/60)41 
but C2 < 3 
| (1. 3x41x Q/60) Cl 
لر‎  —— U SE 
so сўз : کہ‎ al. e (gcc (223-17) 
From equation (2.3-11a), (2.3-12c), and (2.3-17) 
mVx(tanX2 - canto) E .3x41x(£2/60) E Cl of 
ھ۶"‎ c, SM C | С [5м е/суЕ|° 2 
n (A/c2)8 
(2,3213) 
۲ ., _ _AVx(tamk2 رم‎ ۵ 
Stress) = (Sg Eh 8 (2.3-19) 


Freq = ll. pte Say KN (2.3-29) 
A/C 


A e 
1.3x41(£2/60) (stress3) = cl(1 OEE) as 
C7 Freq C2 
C4 = 1.3x41(£2/60) Stress3/Freq (2.3-21) 


! x 


Chord = C = E (2.3-22) 


0۱0 - ct 
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Now that the chord has been determined the bending stress and natural 


bending frequency can be determined from equation (2.3-lla) and (2,3- 
12c) respectively. Having solved for the chord at the base of the blade 


the chord can be determined at mean and tip assuming a linear taper. 


Ab(1 + Areara) 


Abm = eT (2.3-23) 

Assuming constant (A/C?) 
E ue 

Cm = Abm é _ 2: 
TT (2.3-24) 
AO 


At the opposite end from the base 


TE |Arearatio(Ab) š (2.3-25) 
|. (A/C^) 


where С, = Chub (for stator) 


oe Chip (för rotor) 


Now that the chord has been determined at the midchord the aspect ratio 
and pitch can be determined by equation (2.3-26) and (2.3-27). 


Aspect ratio (h/C) = h/Cm 0223-20) 


Pitch = (P/C) Cm (2627) 


With pitch determined the number of blades can be calculated, 


210 


Number of blades (N) = Inteser 
d Pitch (2.3-23) 


The computer program uses the criteria >f having even númber of blades 
for the stator and prime nanbers for the rotor to prevent blade inter- 
action of resonant frequencies. The determination nf whether a number 
is prime is done in subroutine TESTP. Thus havin; solved for the 


number of blades the pitch can be recalculated and the weight 








34 
۳۲۰۲۲ blades. Blades are integrated from the base to end to get 


volume and is multiply by its specfic weight to get the total weight. 
In the computer program this specific weight is used for both the 


blades and disc. Also it is assumed they are solid blades with no 


cooling passages. 


rt 
Weight of blades Ngo J A(r)dr (088095 
r 
h 


x 
- NgbAb Sa - gr) dr 
r 


ya fi دہ‎ ei 
2 


where (X is given by equation (2.3-9) 


Ia this section pitch was determined by a chart of ideal pitch to 
reduce losses. In actual turbine desian the loss data would be 
coupled with stress analysis of che turbine blade attachnent to 
determine allowable pitch. The ty^^ nost commonly used is the fir- 
tree method which is illustrated in figure (2-17). This method of 
selecting turbine blade chord length is just one of many which are 


used in industry. 





Figure 2-17. Firtree 6 
attachments. 
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24 Turbine Disc Desizn 

The method for design of a disc in this section will be based upon 
a constant average stress analysis. The stress relations are those 
proposed in reference (15) and (18). This method will ignore stress 
concentrations and thermal stresses. In actual disc design these cannot 
be ignored since thermal stresses may be a major part of the disc 
stresses and stress concentrations may cause the disk to fail even th- 
ouzn the average stress is well below the ultimate tensile stress of 
the material (UTS). For calculations in the computer program the aver- 


age tensile stress used is that shown in equation (2.4-1) from reference 





m 
۱ vox. (2.421) 
Nb | 2 
No 
Nb | 
: | = Ratio of burst speed to design speed 20022) 
| No 





The value of this ratio is normally between 1.2 - 1.3. The value of 
ultimate tensile stress (UTS) is taken at the desien operating tem- 
perature, higher values producing smaller disc dimensions. This 

value is an input to the computer program. The computer program has 
minimum values for some of the dimensions illustrated in figure (2-13) 
to prevent impossibilities to manufacture and to allow adequate thick- 
ness for attachment of blades 


The following equations are based on figure (2-13) taken from reference 
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Figure 2-13, Turbine disc stress analysis (From Aircraft Turbines 
by J.L.Kerrebrock(15)) 
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for a linear tapered blade, 

For the computer program Wo will be assumed to equal .65 chord 
and fo will equal half of Wo unless otherwise specified. With these 
parameters known Zo can be solved for by equation (2.4-5). With Zo 

Zo Les E^ EI _, GENEE 

Lo م۴‎ 090٥ 2dg0 


known Z(r) can be solved for. 


Zu || bo © 
Lo 
zur) Zo SETS (Es) (2.4-5) 


with a minimum of 0.3 inches, The value of r; iS assumed to equal 
2 ۶ جع وا‎ where rshaft is an input to the computcr program. This 


defines Ti so wi can 1 solve "E in equation (2.4-7). 


UE 
ی‎ CO rt | (2427) 
go 


with a minimum of 1.2 Wo. 





Having determined the dimensions of the disc the weight of the 


disc can be determined by equation (2.4-8) 


Disc weight ٭‎ pf ۵) 0۵9 (2.4-3) 
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Having determined the dimensions of the disc, the burst speed 
can be computed assuming the disc yields prior to bursting. For this 
analysis the disc is divided in half as illustrated in figure (2-19) 
and the stress is assumed constant over area 2A. The force acting 
on this area is composed of the disc load and the blade load, It 
is assumed that the blade load is evenly distributed on the outer 
rim, In actual stress analysis of the disc this assumption is an 
over simplification of the stress, but for preliminary desien 
purposes it will be adequate. A more detailed analysis is covered 
in reference (7). The force due to the disc is now calculated by 
taking the force in the direction perpendicular to the area 2A 
caused by the disc load and summiny it over the disc half and adding 
this to the force caused by tne blade load. This force is divided 
by the area 2A to get the average tensile stress © ta. This stress 
is set equal to the ultimate tensile strength of the material of the 
disc at its operating temperature. The value of Nb is then calculatec. 


by equation (2.4-16) 
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Figure 2-19 
Turbine Disc 
load distribution 
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A = Wi Ti + fauna + ۵٥ (2.5-9) 
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For disc load 


di - Jav (2.4-10) 
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dF = ) كن 12م‎ (5 1 = dor (As ins (2.4-11) 
BEN in sinededr (2.4-12) 
60 haft 
EN ы у? dr 
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wi 





rh 
F = Ta 2 a feta 4 Jana + Wo > Sets | 








| 99 [so Tshaft 
、 3 3 
_ (2N 2 Wi(r; rar) ۲ „ch are meter _ To 3) 
60 / نہ‎ 3 
ы ES 
2, Re q 
رما‎ e - 22 ^ e 
0 e r dr is integrated numerically 4 =. E EY 
1 ave 
E 
Fd; 
ete (due to disc) = - (ONES 
For 5lade load 
N 
Sta(due to blades) = шы) cb given by equation(2.4-4) 
(2.4-14) 
Average disc stress 
CSta(Disc + blades) = Pdisc * Fblades (2.4-15) 


2A 


Set cita = ULT (ultimate tensile strength of material) 


The burst speed Nb can now be calculated. This value will 
Slightly higher than that in dquation (2.4-2) due to approximetions and 


arbitrary values selected for some of the disc dimensions. 


12 


- 83 e —.. 4 
rh A 2 ab 
1 ۱ 


= 
ae 2 SE 











(2.4-16) 





41 

2.2 Having determined the blade and disc dimensions the axial 
۶۳۲۲۱ Of the turbine can be calculated. To do this a criteria for 
Spacing between blades must be selected. From reference (6) the 
values of spacing between blades is 0.2 to 0.5 of the axial chord 
length. The longer dimension will tend to reduce the flare which 
would prevent the gas flow to separate from the turbine wall. It 
would also decrease the chance of inducing vibrational stresses. 
The shorter dimension would decrease the size and weisht of the 
turbine. For the computer program a value of À the average axial 
chord at mid-heisht was selected. The computer program will then 
calculate the axial lensth and weisht of the shaft from station 1 


of the lst stage to station 3 of the last stage. A value of fo 


was selected as the stagger angle in the computer program, 


Figure 2-20 Axial 
Spacing between blades 
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Figure 2-21 Stagger Angle 
of Turbine blades 
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Chapter 3 42 
Performance Estimation 


21 The design and off-design characteristics are very important 
in a complete analysis of a gas turbine since most turbines do not 
Operate only at their design point. Also cycle calculations without 
actual loss data can not accurately determine the efficiency of a 
turbine, so this analysis can give a more accurate assessment of the 
efficiency that can be achieved in the turbine. 

The Ainley-Mathieson Method for performance estimation was 
selected over other methods because it seems to have a wider applic- 
ation and can be easily used for off-design performance calculations. 
The modifications proposed in reference (8) have been applied. These 
modifications can reduce the error from *3;4 to $27 for most turbines 
and are based on modern turbines while the data by Ainley and Math- 
ieson were accumulated over 20 years ago on turbines that existed at 
that time. The calculation procedure for off-design performance is 
that proposed by Horlock in reference (12). A description of the 
Improved Ainley Mathieson Performance Estimation Method will be pre- 
sented first followed by a description of off-design calculations which 
use loss coefficients to determine gas properties, isentropic efficiency 
and static efficiency at design and off-design points. 

The computer program used pressure loss coeíficient Y = (Po:r- Poz) 
/ (Poz -P; ) for the performance analysis. This cocíficient is a iunc- 
tion of aspect ratio (H/C), pitch to chord ratio (P/C), velocity, Mach 
number, trailing edve thickness, thickness to chord ratio (t/C), 


Incidence angle, tip clearance, Reynolds number, etc. This coefficient 
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Loss Coefficient 
„ Poı -Po2 
Po2-P2 
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Figure 3-1. Analysis of losses in flow 
through a row of turbine blades( From 
D.5.Ainley and 3.C.R.Mathieson(1)) 
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is broken down into three parts, Yp (profile loss), Ys (secondary loss), 


and Yk (clearance loss). Figure (3-1) from reference (1) shows the 
first two, Yp and Ys in which the annulus loss is grouped into second- 
ary losses. The modification of the Ainley-Mathieson method mostly 
affects the secondary and tip clearance loss coefficients, but has an 
affect on profile loss Yp for Mach numbers greater than 1, and a Rey- 
nolds number correction to the profile and secondary loss coefficients. 
Do? The pressure loss coefficient Y is calculated in subroutine 
PLOSS. Figure (3-2) shows the sign convention used in this subroutine, 
where 4 angles are blade angles and OX angles are gas angles. The ang- 
les as shown are in the positive direction. The (') implies relative 
angles. This sign convention and notation applies only to subroutine 
PLOSS in order to have a standard form for stator and rotor. 

Figure (3-3a) and 3-3b) coupled with equation (3.2-1) is used 
to calculate the profile loss coefficient Yp for O incidence, where 
£i Q values come from figure (3-3a) and Bur = (<> comes from figure 
(3-3b) The pitch to chord ratio (P/C) is the argument for entering 
figures (3-3a) and(3-3b). The correction for thickness to chord 
ratio (t/C) is done in equation (3.2-1) along with the correction for 
the ratio of inlet blade angle to exit gas angle 1/09). The value 
obtained is Yp for zero incidence (i0). 

For off-incidence calculations figures (3-4) a,b,c, and figure 
(3-5) are used. Figure (3-4) a,b,c determine the stalling incidence 


so that figure (3-5) can then be used to correct Yp for off-incidence 


Cases, 
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2 
Yp = ar + (8,/ 04) [Yp(4,- 9 07 0) 
155) (B,/ 06) 
0.2 


Xa 
Of P/ C». 75) = (al Col P/ C=. 75)) (3.2-2) 


Figure (3-4b) is entered with the argument of pitch to chord 
ratio (P/C) . The ratio of exit gas anale OX, to exit 2091606 { P/C=.75) 
is then determined. It: is.now possible: to determine 0£(P/C+. 75), Having 
determined(((P/C=.75), figure, (3-409 T new entered with do (P/C=. 75) 
to get the stalling incidence ig (P/C=.75). Figure (3248) is now entered 
vith exit gas angle (and P/C to get A ig = is - is (P/C=.75). ig can 
then be: solved for. ilaving incidence angle i and is the ratio 
(i/ig) for entering figure (3-5) to determine ratio Yp/Yp(i=0)4s known. 
Having previously solved for Yp(i#0) in equation (3.2-1) Yp is now 


known for off incidence. 
2 
Ys = .0334(C/H)( cosd}/ cos) Z (3.2-3) 
For secondary losses Ys is computed by equation (3,2-3) 


from reference (8). Where Z is the Ainley loading parameter from 


reference (2). 


C cos + 2 
ai cos $ ue. кош 
Km = (W- X9/2 (for stator) (3.2-5a) 


(m = IW EL (for rotor) (3.2-5b) 


-> 





47 





-20 E ^ V 
Values of (Xy 60 









EN = iq 
i, (P/ C=. 7590 


| | | 
0.4 05 0.6 0.7 0.8 0.9 ] O 
Pitch/Chord (P/C) 


| Stalling, 
Incidence 
is when 
P/C = 0.75 
20 


LU 


۰ mm E eege, 


y 0S сол 0.0 0.4 0.3: 1.2 
(c) 


Figure 3-4. Of f-design performance of cascades of turbine 
blades. Re = 2x10“: °` .€a)variation of stalling incidence 
with nitch-chord ratio;(b)variation ofX,with p/c;(c)stall- 

| inc incidence of turbine blade sections when p/c = 0.75.(. 
From D.3.Ainley and G.C.R.Mathieson(2)) 
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Figure 3-5, Variation of relative profile loss with 
relative incidence. (From D.G.Ainley and 5.C.R. 


Mathieson(2)) 
Li = 2(tan04 + tan (7) (3.2-ба) 
ч / 
E = 2(tanQGs tang (3.2-6b) 


For tip clearance losses Yk is computed by equation (3.2-7) 


from reference (3). 


3 
Z 


Yk = لت‎ (E) 
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For equation (3.2-7) 


d =» 


k = 


C = 
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Fizure 3-0. effect of trailing edge thickness 
on blade coefficients. (Fron D.3.Ainley and 
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deternined Yp and Ys they can be corrected for Mach 


number and aeynold's number from equation (3.2-3) and (3.2-9) 


respectively. 
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(Yp) corrected = Yp(1+60(Mn-1)*) (3 228) 


for Mn Pl 


Re 7°? 


23 (3,229) 


(Yp + Ys) corrected = (Yp + E 


Now the total loss coefficient Yt can he computed. by equation 
(3.2-10). This coefficient is now corrected for thickness of trailing 
edse to pitch ratio (te/P) from figure (3-6) to set the total loss 
coefficient for the blade row. Illustrations of an example calcula- 


tion can be found in reference (2) and (12). 


Yt = Yp + Ys + Yk (3.2-10) 


IW 
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B.) In this section it will be explained how the loss coefficient 


Yt solved for in section 3.2 is coupled with compressible flow tables 
to estimate performance. fhe method explained is that outlined in 
reference (12). Another method is cxplained in reference (2). Cal- 
culations will be carried out only at the mean radius which should 
give a fairly accurate assessment of .the turbine performance. 

[he calculations are started at the inlet of the turbine and 
worked through each stage where the exit conditions at one blade 
row is the inlet to the next blade row until the turbine exit of the 
last stage is reached. 

For inlet conditions the values for Toi and Poi will remain 
constant, the mass flow and corrected speed 0 [fo3/ ۷/۳۵: مه(‎ will 
will be varied. This will allow the values of efficiency verses 
pressure ratio and уне verses pressure ratio to be cal- 

Poi Poi /dp 
culated as illustrated in Figure (3-7). 
3. 3a Flow through the nozzle. 

(1) Select the mass flow. 

(2) Assume an exit Mach number. 

(3) With known blade angles subroutine SIOCRA is entered as 
in section (2.2) only this time the blade angles will be used to det- 
ermine the exit gas angles. 

(4) From compressible flow tables for the assumed Mach num- 


ber the values of (Y./Po2) and (v / 1020 are determined, which will 


then determine V 9 and و1‎ 


Vo = [Toz (v,/ [To2) (3.3-1) 


SE 
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Static efficiency vs Pressure ratio 
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Figure 3-7, Performance Curves for Axial Turbine (a) 

Issentropic efficiency vs pressure ratio. (b) Static 

efficiency vs pressure ratio. (c) Corrected mass flow 
vs pressure ratio, 
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Subroutine PLOSS is now entered to set Yt (pressure‏ )5( 
loss coeffiecient).‏ 
Po2 can be determined from known values of Yt and‏ )6( 
(Po /Poz) from equation (3.3-3)‏ 
Poz = Po1/(Yp(1-P7/Poz) + 1) (3.3-3)‏ 


(7) The mass flow parameter Q2 can now be determined. 


my To 
E A» Po ( ) 
where À> = و۸۰‎ 5 (3.3-5) 


(8) From the compressible flow tables for the value of Q, 
the Mach number can be determined, and if this value is close to the 
assumed Mach number the inlet conditions to the rotor can be solved 
for the new values of V یرہ‎ Mach number, and T will be used 
to iterate the process until the assumed Mach number approximately 
equals the calculated value. Also as the Mach number approaches (1) 
a calculation will be performed to see if the row is choked. (ie; 
maximum flow for given parameters). This will be covered in section 
EC. 

Now that the exit Mach number has been determined the 
exit conditions of the nozzle can be determined, which will be the 
inlet conditions to the rotor. Where ۷ is determined from equa- 
tion (3. 3-1), T, from equation (3.3-2) and P) from compressible 
flow data. 

B ^ Po2(P5/Po2) (3.3-6) 

8 p (3.36a) 

RT, 
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3000 Flow through the rotor 


(1) ۷ 4 V „known the relative rotor angle 4, can 


be calculated. 


Um» 
ta 3 = tan 2 (3.3-7) 
$ x E و۷‎ ء٤4‎ 


where Um = tm GJ 





t 
and determine the incidence anale (i,) where (5, is the blade angle 


and i45, = => 


(2) The relative velocity W and relative total temper- 


ature, can then be calculated. 


m 
Wo- (Vacos 2 J) (F taf 42) (3.3-8) 
2 
Tozrel= 了 + W2 (3.3-9) 
230JCp 


(3) W2/4 To? 6} can be calculated, which will be used to 
enter the gas tables to get P»/Po2 „eı which will allow De em 
to be calculated. 

(4) The same procedure of assuming exit Mach number and 
reiterating until the assumed Mach number equals the calculated one 
is followed. :1م503‎ 5 5 W4,M3rel,and Datei will then be known. From 
this relative Mach number (Py/Po2r.]) can be determined, which will 
enable P4 to be solved for. The absolute parameters (Х„, Va 
Po3, T4 and To3 can then be calculated. 


Um 
t = t + = 
an (X4 an, 1,002, (3.3-10) 


where Um = Can لب‎ 


١ 4 


- 


mt 





(W_cos (| ta £43) (3.3 jus 
3 ¿Cos f3 + n ox. کا‎ 


<= 
l 





I4 = Toa D 00 
3 nes 230JCp ` 
Poa = P3(To3/T3) Y-T (3.3-13) 


(5) If this is the last stage, efficiency Z can be calculated 
by equations (3.3-14) and (3.3-15), If not the exit conditions will 
determine the inlet conditions to the nozzle of the next stage. This 
process will be iterated until the parameters have been determined for 


every stage in the turbine.. 


_ (Toi - Toe) oo 
Wee 5 iin - ix d сы) 
Hes = ا سم لظف‎ (3.3-15) 


roi |1 - (Pe/Poi) ¥ 


Boc Determination of whether a row of turbine blades have choked 
is carried out when the Mach number approaches (1). For the computer 
program this process will be carried out when the Mach number is creat- 
er than 0.9. 

(1) Select a slightly hizher value of Q m... calling ict 
ON 

(2) With this Q* enter the compressible flow tables to get 
corresponding values of P/Po and solve for Pojn /Poout assuming Yt 


remains constant. 


= 





20 
(3) Now calculate the new mass flow m, Tf this is less than m 


the flow has choked. Yt and Q will be assumed to remain constant with 
the flow going supersonic at the exit of the blade row. If m' is not 
less than m calculations will be resumed as was illustrated in section 
Misa) and (3.3b). 

(4) Enter compressible flow tables with Q for supersonic flow 
to get corresponding values of (Pout/POoout)»s Vout / [Toout aE MES 
From equation (3.3-3) determine Pout. With equation (3.3-4) Aout can 


be determined. Since Aout = A cos(X , (can be solved for in 


annulus 
equation (3.3-16). The calculations are then resumed in the same 
manner as in section (3.3a) and (3.3b). The validity of this method 


depends on whether Q and Yt remain constant after the turbine blade 


row has cnoked. 


(X= cos mToout —_ (3.3-16) 


AannülusPoQ 


b 


i? 


RT 


bd 
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Chapter 4 


Operation of Computer Program 


Chapters 1, 2, and 3 described how the calculations are‏ ات 
performed in the computer program. This chapter will show how the‏ 
flow of information is passed throuzh the computer program to get‏ 
the desired results. This information is not needed to use the‏ 
computer program since that is covered in Appendix (2).‏ 

4.2 The steps of calculations in the computer are in the same 
sequential order as illustrated in chapters l through 3. This is 
shown in the simplified block diagram figure (4-1) 

4.3 The method used in the computer program to interpolate 
curves is a 4 point Lagrange Polynomials with a worse case of a 
linear interpolation due to break down of Lagrange Polynomials for 
monatomic increasing or decreasing function. Subroutine FIG and 
BK are used to do this interpolating. Subroutine FIG selects the 
4 data points closest to the argument X for passing to subroutine 
BK which actually does the interpolating to get the ordinate Y as 
an Output. If the argument exceeds the limits of the data points 
the value closest to the argument X is used to get the ordinate Y. 
The reason this method was chosen over finding formulas by least 
squares or simular methods is this method can be adapted to chan- 
ges in data without solving for new formulas or coefficients. The 
Lagrange Polynomials uses a monatomicly increasing base to get the 


Ordinate Y. If the base is not monatomicly increasing such as 
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values of Mach Number vs Ps/Po, subroutine FIG rearranzes the base so 
that the base is monatomicly increasin» and the ordinate values are 
in corresponding relationship with the base. For values of base such 
as (m]Toi)/APoi ) which both increase and decrease the data »oints are 
divided into two parts with the part that increases as one set of data 
points and the part that decreases as another. [he ordinate does not 
have to be a monatomic function. The problem of increasing and 
decreasing bases can be enccuateéred in the performance estimations be- 
cause the base and ordinate are interchanged in some of the steps. 
Figure 4-2 illustrates how subroutine FIG and 3K are used to inter- 
polate data. 

4,4 Figure (4-3) illustrates the method the computer program uses 
to determine the flow area at station l. In most cases for a value of 
epsilon (€) = .001 the convergence is reached in less than 10 iter- 
ations. This method is used for solving for the exit flow area and 
also for station 2. The method for solving the other parameters 

were will documented in chapers 1 thru 3 and will not be repeated 
Шеге. 

225 From output generated by the compüter program, it is felt that 
the computer program can be very useful from early preliminary desien 
to start of detailed design, since it has the versatility to accept 
more input data as the designer firmly sets the characteristics of the 
Turbine. As for cooling requirements it is felt that compensations 


in the input data could account for losses attributed to cooling. 
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BLOCK DIAGRAM OF SUBROUTINE FIG AND BK FOR INTERPOLATING 90 
FUNCTIONS f(x) from MAIN PROGRAM OR OTHER SUBROUTINES 
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Axial Velocity Vx at Station l1 For Constant 


| 
Hub Turbine 


Figure 4-3. Flow diagram 
for determining flow 
area at station l 
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APPENDIX I 


Symools and Notation used in this Text. 
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ti 





ho 


H/C 


Po 


P/ C 


To 


SYMBOLS USED IN TEXT 


Area 
Chord of turbine blade 
Specific heat (btu/lbm?R) 
Constant in Newtons Law 32.174 (ftlbm)/(ibfsec?) 
Blade height 
Stagnation Enthalpy 
Aspect ratio 
Incidence angle 
Moment of inertia 
Energy conversion factor 778.16 (ftlbf)/(btu) 
radius of gyration 
Mach Number 
vx/(_2. goRTo) * 
+1 
Mass Flow 
Number of Stages 
Number of Staze 
Stagnation Pressure 
Static Pressure 
Pitch to Chord ratio 
Gas constant (£t1b£)/ (Lbm*a) cp 
Corrected mass flow mjTo/PoA 
Stagnation Temperature 


static Temperature 
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Ti 


SM 


Vx 


Wo 


Wi 


EC 


N 


Zi 


Zo 


Z 


on 
۳۹ 
мы 
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aR eta ® A 


Thickness of Outer Rim 

Thickness of Inner Rim 

Non Dimensional Section Modulus 

3lade Speed 

Velocity 

Axial Velocity 

Relative Velocity 

Width of Outer Disc Rim 

Width of Inner Disc Rim 

Pressure Loss 

Radius 

Radius 

Thickness of Turbine Blade 

Thickness to Chord Ratio 

Section Modulus 

Thickness of Mid Section of Turbine Disc at Inner Radius 
Thickness of Mid Section of Turbine Disc at Outer Radius 
Thickness of Disc as a Function of Radius 
Absolute Gas Angles (except section 3.2) 
Relative Gas Angles (except section 3.2) 
Stagger Angle 

Ratio of Specific Heats Gamma 
Efficiciency 

Stress 

Density 


Psi Loading Coefficient = 0ع‎ 


Uh 


0 


bh 


(1 


—— 





(R Reaction 
EP speed (rpm) 


C4) Speed (rad/sec) 


Subscripts 

i Inlet 

> Exit 

t Total 

S Stage 

p Poiytropic 
S Static 

t Tip 

m Mean 

h Hub 


C constant 
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APPENDIX LI 


User's Guide To the Computer Program with Comparison of 
Turbine Desizned with this Computer Program to 
one which was. Designed by Airesearch. Sample 
TR and Output of Turbine 
Designed by Computer 


Program 
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11-1 The computer program can be used to calculate p..ysical and 
thermodynamic properties of an axial flow turbine.. It nas an optional 
ability to do performance calculations. The optional phase is very 
expensive to operate since it does both design and off design calcu- 
tations. From Operating the computer program and comparing the cycle 
calculations with performance calculations the cycle calculations 

are within 3% of performance calculatians at cesign point. For the 
early preliminary phases of turbine desien it is recommended that 
just the cycle calculations and mechanical design phases be carried 
out with the computer program punching output cards to have data 

to perform the performance calaulations at a later time 

Á.second prosram which reads in the punched output cards in the 
same order as they were punched in the first program is provided to do 
only performance calculations. This allows the performance estimations 
to be carried out with runniding the complete program again. 

For future references program #1 will refer to the complete 
computer program while program #2 will refer to the program which only 
does performance estimations. The complete computer program can not 
be operated in computers which have less than 210 k of memory. 

11-2 Pages 7] to 76 show how the computer cards are to be punched. 
The lst card is a label which will be printed out in the same content 
as it was entered on the computer card.. Cards 2 thru 4 are the input 
data cards of known turbine parameters. Cards 5 thru NNN are the 
compressible flow tables and viscosity of the gases with correspondins 
temperatures. Starting from card 5, lst the Mach Numbers are read in 


» 2nd the corrected mass flow m/To/PoA, 3rd the préEsure ratio P/Po, 
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70 
۱۱۳۱ ۱6 corrected velocity v/{ To, 5th the temperatures which correspond 
to the viscosity data, 6th the viscosity data. Pase // is a set of 
compressible flow tables and viscosity data for air. 
II-3 Pages 73 thru 82 are comparison data of computer program designed 
turbine with one which was designed by Airesearch reference (20). AS can 
be seen from the data some of the parameters are different because of 
lack of knowledge of some of the design parameters which was used in the 
design of the Airesearch Turbine, also some simplifications are 
incorporated into the computer program for solutions of output para- 
meters. 
11-4 A list of sample output and input is provided from page 83 to 
93 for turbine which was compared with Airesearch turbine. As can 
be seen on page 91 of output data, the performance data is very sim- 
ilar to the cycle calculations at design point. This justifies the 
argument of doing only cycle calculations and mechanical calculations 
when off design data is not required. Figure II-3 is a plot of the 
data pcints from performance calculations. Where some of the curves 
are discontinuous is caused by a limit of 10 iterations for solution 


at a blade row.. 


i 





Input Data Cards 


* (Must be specificd) 


71 


# (If value of parameter is not known place value to right of 1 


on computer card) 


Bard 1 


Column 1-80 


Alphanumeric label (ie, "Gas Tur- 


bine number 1") 


# Blank Card 


Card 2 
| Column 1* Ml (Integer) 
Q compressible flow data to be 
entered, 
1 No compressible flow data to be 
entered. 
| Column 2 * M2 (Integer). 
x O Power in horse power specified. 
1 Pressure ratio specified. 
Column 3 and á Blank 
Column 5 * M (Integer) Type of turbine. 
l Constant hub 
2 Constant mean 
3 Constant tip 


Note 


Where instructions inform user to place dummy values on 


cards these values are used 


to inform computer that they 


are not specified and to perform optional calculations to 


determine the parameter. 
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Column 6-10 


(leave blank if Ml = 1) 


Column 11-195 


(leave blank if Ml = 1) 


Column 16-20 * 


Column 21-25 * 


Column 26-35 


(lbm/sec ft) 


Column 36-45 * 


Column 46-50 * 


Column 51-55 * 


Column 56-60 * 


72 
KLM (Integer right justified) 
Number of compressible flow el- 
ements to be entered. (ie num- 
number of Mach Numbers to be en- 
entered) ( Maximum 100) 
KIM (Integer right justified) 
Number of Temperatures at which 
values of (viscosity) x 10 No 


sel ft. 
is defined. 


Gamma X (FS 3) 

CP (3tu/lbm PR) 

(F5.3) 

Power (Horse Power) ( F10.1) if 
M2 = 0 this must be specified 

# 9000. 

W (mass flow) (lbm/sec) (F 10.3) 
TC (thickness to chord ratio) (F5.3), 
normally about # 0.2 

TET (thickness of trailing edge to 
pitch ratio) (F5.3) normally about 
# 0.02 

TIPCLA (Tip clearance) (in) (F5.3) 


normally about # 0.015 
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Card 3 


Column 1-5 * 


Column 6-10 * 


Column 11-15 * 


Column 16-20 


Column 21-30 * 


Column 31-32 


Column 33-39 


Column 40 * 


Column 41-45 * 


73 


p PSI (F5: 3) 


Normally between 1.5 and 2.9 

# 2.0 

RTE (Max Tip raduis in ft at exit) 
(E533) 

Ratio T (Ree/Rne) Ratio of Tip to 
hub at exit (F5.3) # 1.5 

(X hub (exit) Reaction at hub (F5.3) 
normally between .00 to 0.4 # .05 
OMEGA (RPM) (F 10.0) 

KLMi (INTEGER) Right justified if 


Compressible flow data provided 


(MA Toe 
tor A 


is a maximum in (for 


the value at which 


sample problem KLMl = 50) (Maximum 70) 
Blank 

M4 (Integer) 

O performance calculations to be per- 
formed. 

1 No Performance Calculations to be 
performed. 

SPECWT (Specific Weight of turbine 
disc and blades (1b/in2) (F5.3) 


# normally about 0.3 


t 


> 
- 





Column 46-50 * 


Column 61-70 * 


Column 71-80 * 


Card 4 


Column 1-5 * 


Column 6-15 


Column 16-20 


Column 21-30 * 


Column 31-35 


Column 36-45 


74 
Y MOD (Young's Modulus psi) (f 15.0) 


Normally about # 30000000. 

Cl (Value of allowable blade stress 
in psi) (F 10.0) Normally about 

# 50000. 

C2 (Value of allowable blade steady 
stress in PSI) (F 10.0) Normally 


about # 100000. 


ENT (Value of total efficiency 
desired. # .9 
Radius of Shaft in ft (F 10.5) # - 0.1 
if not known 
P ratio (Pressure Ratio) (F5.3) 
if M2 = 1 this must be specified # 1.5 
ULTSRE (Ultimate Tensile Strength of 
Material of Disc (psi) Value for disc 
stress is calculated from this. (F10.0) 
A ave = .75 ULTSRE 
2 
012) 
Normally about # 130000. 
ALPHAH (Exit angle in radians) (F5.3) 
Normally between (0.0 to - .3) #-.1 


EMACH = Vx Normally Between 





2 
soRTo 
0 +1 


2-10 #-. 5 
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Column 46 


Column 47-48 * 


Column 49-50 


Column 51-55 * 


Column 56-60 


M3 (Integer) 75 
O Normal Blades 


1 Shrouded Rotor blades. 0 


NKL (Number of turbines to be de- 
signed by computer program. (In- 
teger right justified) All sets of 
cards except for last turbine des- 
ign will have O in column 48. The 
last set of cards will have the núm- 
ber of sets that were entered. This 
allows the computer to stop after com- 
pleting the last turbine design. 
N (Number of stage Integer right jus- 
tified) (Maximum 10) If this value 
is not known, enter O and computer 
will solve for number of stages. ¿0 
AREARA ( F5.4) 
Tip to base area ratio of turbine 
blade normally between „25 to 5 

# .25 
TOO (Gives radial thickness of outer 
rim of disc. by To--=Wo/TOO normally 


between 2 to 4 (F5.3) 4 - 5 


76 
Cards 3-NNN 


These cards are optional. If performance calculations are 
to be performed they must be furnished. If more than one turbine 
is to be designed compressible flow data must only be furnished 
for 1 st turbine and the value of Ml on card 2 for each following 
turbines will be "1", If the compressible flow parameters are dif- 
ferent a set must be furnished. A set of compressible flow data and 
viscosity data for air are provided with the computer program. 
Card 5 - N RMAC (Mach Numbers) (20F4.2) 

same number as KLM 


° O 
Card (N + 1)-K WTAP (10F8.5) ANTO 1o R. 


PoA sec 1bf 
Same numbers as KLM 


Card (K+1) -J PSPT (Ps/Po) (10F8.5) same 
number as KLM 

Card (J+1)-L VELTOT (V/ J To) (10F8.5) (ft/sec) 
same number as KLM ١ 

BASAL + 1) - LL TEM (°R) (20F5.0) same number as 
KTM temperature values that cor- 
respond with viscosity data 


Card (LL + 1)-NNN V (Viscosity of fluid (lbm.sec/ft)x 


1000) CORSO) 


=, 


CS 








Compress 


< —— E  . — ” نب‎ 


- 
MACH NO w* SORT (T) / (AP) 


...0.02000 
0.04000 

. 0.06000 
0.08000 
0.10000 
0.12900 
_— _0.14000 
0.16000 
0.17000 
0.20000 
0.22000 
0.24000 
۔۔‎ 0 
0.28000 
--.0.30000 
0.32000 
0.34C00 
0.36000 
0.38000 
0. 40600 
___0.42000 
0.44000 
0 0 
0.48000 
0.50000 
0.52000 
0.54000 
0.56000 
0.58C00 
0.60000 
_____ _0 0 
0.64000 
_ 0.466000 
0.68000 
0.70900 
0.72000 
0.74000 
0.76000 
0.78600 
0.80000 
Lo... 0.482000 
0.84000 
LL. 0.486900 
0.88000 
. 0. 90000 
0.94000 
0.96000 
—..... 0.98000 
1.00000 
1.02000 
1.04000 
1.06030 
1.08000 
1.10000 
1.12000 
1.14000 
1.16000 
- .. 1.18000 
1.20000 
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0.01811 
0.05296 
0.05388 
0.C7171 
0.08945 
0.10706 
0.12448 


0.14154 


0.15885 
0.17562 
0.19226 
0.20849 
0.22446 
0.24049 
0.25612 
0.27135 
0.28590 
0.30026 


. 5 


0.32763 
0.34118 
0۰ 35 3 6 
0.36627 
0.37821 


0.38948 | 


0.40061 
0.41108 
0.42114 
0.43067 
0.43979 
0.44838 
0.45669 
0.46421 
0.47159 
0.47796 
0.48433 
0.49010 
0.49560 
0.50039 
0.50459 
0.50815 
0.51170 
0.51479 
0.51734 
0.51948 
0.52123 
0.52257 
0.52337 
0.52391 
0.52418 
0.52391 
0.52337 
0.52257 
0.52135 
0.51989 
0.51801 
0.51599 
0.51358 
0.51076 
0.50795 
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COMPARISON OF COMPUTER PROGRAM TURBINE DESIGN WITH 78 


AIRESUARCH 9000 HP FREE VORTEX AXIAL FLOW TURBINE 


Inlet Conditions 
Sas Air 
Poi(psi) 312.5 
Toi( OR) 1767.3 
m(1b/sec) 119.421 
Design Parameters 


bois 


Xe 0.0 
Assumed Design Parameters 
Mer 0.9163 
G»p(btu/1b?R) .2715 
Gamma X 1723353 
Results 
Figure II-l is a comparison of the turbine wheel dimensions, 
Figure II-2 is a comparison of design point performance calcu- 


lations of total and static efficiencies as a function of tip clearance. 
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COMPARISON OF OUTPUT DATA FROM COMPUTER PROGRAM WITH 


Parameter 
Toe (FR) 


Mori 


Poe (psi) 


hub(station 1) 
radius(ft) 
eos? 
QC, (des) 

tip(station 1) 


radius(ft) 
(3, (des) 
A; (des) 


Hub(station 3) 
radius(ft) 
dee 
A , (des) 
OC, (deg) 


Tip(station 3) 


radius(ft) 


O, (dez) 
OC, (dea) 


d cf (Psi) 


wheel Weisht(lbf) 


at exit 


AIRZSEARCH GAS TURBINE 


Airesearch 


1+2 


136.0 


Velocity Diasrams 


4333 
-62.02 


0.0 


0ء 
67.15- 


0.0 


.4383 
-60.95 


0.0 


6215 
-68.26 
0.0 
37,400. 


13.7 


Computer Prosram 


1554.0 


.429 
-63.96 


059 


529 
-70.93 


0.0 


22 


-63.96 


.036 
< 
0.0 
37,396. 


21.46 


| " 
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Fonelusions: 

Most of the descrepancies in parameter values can be attributed 
to different input values for the computer model. This is caused by 
not knowing the assumptions or definitions that were used in designins 
the Airesearch Gas Turbine. From this limited comparison it is felt 
that the computer model can be used for designing marine gas turbines 
Since it has the versatility to change parameters which are attributed 
to specific turbines. For high performance turbines (aircraft) it - 
should be compared with known aircraft turbines to verify the validity 
of the computer program and adjust the variable parameters to match 


them. 
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Figure II-l. Comparison of Computer کچ ہے‎ Turbine disc 
with Airesear Designed flosed Cycle 9000 HP Free 
Voftex Turbine Disc. 
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TOTAL AND STATIC EFFICIENCIES 
95 . A VERSES 


7 TIP CLEARANCE 
£ £ 


70 


(total) 
94 Computer design His 


Airesearch Turbine (static) 


89 


87 
86 


تک 





© 0/6 0.30 
TIP CLEARANCE (IN) 


Figure II-2, Comparison of Efficiencies verses 
Tip Clearance of Performance Calculations. 
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